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Overview 


A bright impact flash detected by the NASA Lunar Impact 
Monitoring Program in March 2013 brought into focus the 
importance of determining the impact flash location. 

A process for locating the impact flash, and presumably its 
associated crater, was developed using commercially 
available software tools. 

The process was successfully applied to the March 2013 
impact flash and put into production on an additional 300 
impact flashes. 

The goal today: provide a description of the geolocation 
technique developed. 
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Impetus for locating flashes 


• Bright impact flash on 17 March 2013 

- Thought to have produced a fresh crater detectable by 
Lunar Reconnaissance Orbiter (LRO) 

- LRO needed an accurate search area to begin looking for 
the crater 

• Geolocation workflow was developed to provide 
coordinates to LRO 

— Rough process used for 17 March flash LRO found the 
crater! 

— Refined process put into production on 300 impacts 
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Linking flash to fresh crater 


• Idea of linking observation of an impact flash with its 
crater is an appealing one 


• Provides “sanity checks” for 

- NASA photometric calculations 

— Crater scaling laws developed from hypervelocity gun 
testing 

• Luminous efficiency estimates can be made by 
combining flash and crater measurements, assuming a 
crater scaling law 
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Terminology 


Geolocation 


Georeferencing 


Process of identifying the real- 
world spatial location of an object. 


Method used to associate an image 
with a map of real-world locations 
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Geolocation workflow 





Jacobs, ESSSA Group/MEO/D.E. Moser 


Lunar Impact Workshop, 2-3 June 2015 


6 











JACOBS 

ESSSA Group 


1. Create video segment 


• Identify video frame(s) 
containing impact flash 

• Use YirtualDub to create 
short video segment 
containing impact flash 


VirtualDub 1.6.19 - [4jan2008_SDJmpact2_llh42m38s.avi] 
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2. Determine flash centroid 



• Import video frame into MaximDL 

• Calculate flash centroid location in 
image coordinates, (x,y) in pixels 

• Export video frame as TIFF image 

Record {xjiyj) 
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ArcMap 



• Part of the ArcGIS 
software suite 
developed by Esri 

• Off-the-shelf solution 
for georeferencing 
lunar impact imagery 


Q Lunar_LRO_basemap_30Mar2009.mxd - ArcMap 



ArcMap interface with lunar basemap 
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3. Define & setup basemap 


• Load lunar basemap into ArcMap 

- Source: http://wms.lroc.asu.edu/lroc/view_rdrAVAC_GLOBAL 

- LRO-created orthographic projection of the lunar surface 

- resolution = 32 pixels/deg, center = (0° N, 0° E) 


• Adjust basemap display to account for lunar libration 

- JPL Horizons sub-observer lat/lon used to adjust map center coordinates 


LRO basemap 
centered at 
(0°N,0°E) 



LRO basemap 
centered at 
(7.2°N,7.0°W) 
to georef flash 
2006 Sep 28 
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4. Georeference flash image 



General steps for georeferencing an image 

1 . Overlay image on the basemap 

2. Link known positions in the image (x,y) to 
known positions in the map (x \y ’) using control 
points 

3. Save the transformation used to align the image 
to the map 

4. Save fit error estimate for use in uncertainty 
determination 
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4. Georeference flash image 



Start with video frame 
TIFF 
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4. Georeference flash image 




Image's brightness & 
contrast adjusted to 
emphasize prominent 
features 


Overlay video frame 
TIFF on basemap with 
libration adjustment 
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4. Georeference flash image 




1 control point added, 
linking image to map, 
using ArcMap's 
"Georeferencing 
Toolbar" 
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4. Georeference flash image 



2 control points 
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4. Georeference flash image 



3 control points 
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4. Georeference flash image 




Typical control points 

• High albedo craters 

• Transitions between mare & highland 


20 control points 
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4. Georeference flash image 



Final georeferenced 
impact image 


Jacobs, ESSSA Group/MEO/D.E. Moser 


Lunar Impact Workshop, 2-3 June 2015 


18 


JACOBS 

ESSSA Group 


4. Georeference flash image 



• ArcMap uses control points and a least-squares fitting 
algorithm to determine a 1 st order polynomial 
transformation 


x' = Ax + By + C 

y’ = Dx + Ey + F 

• Transforms image coordinates (x,y) in pixels to basemap 
coordinates (x y ’) in meters 

- (x’,y ’) is the orthographic projection of the 3D Moon onto a 2D 
plane 

• Parameters A-F are determined by the control points: they 
scale (m x , m ), shear/skew (k), rotate (t), and translate all 
coordinates in the image to map coordinates 

A = m x cos t B = m y (k cos t — sin t) C = translation in x direction 

D = m x sint E = m y (k sin t + cost) F = translation in y direction 

Lunar Impact Workshop, 2-3 June 2015 
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4. Georeference flash image 



• Applying the transformation to each control point 
yields a residual error s 

• ArcMap displays s for each control point and 
calculates the root means square (RMS) error 

V n c2 

RMS error = 1 n = # of control points 

-\j n 

• The RMS error is saved for subsequent uncertainty 
calculations 
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5. Transform flash coordinates 



• ArcMap transforms image coordinates to map coordinates, 

(x, y ) (x \y ’) 

• The same transformation is used for the flash centroid, 

(Xf,y f ) (x/,y/) 
determined in Maxim DL 


• Custom code read in the world file containing the 
transformation parameters. Output is flash location in 
meters on the orthographic projection plane. 
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£sss°i? 6. Determine flash location 


• Input flash location to 

ArcMap’s “Go to XY” tool 

• Read & record selenographic 
coordinates (A,, cp) transformed 
by ArcMap 

• Place marker at flash location, 
add point to database and 
shapefile 
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Geolocation workflow 
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